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A method for calculating the Si0 2 diffusion coefficient is proposed and a particular example is described using 
a generalized kinetic dependence of the process of dissolution of silica grains in melted glass. 


The rate of glass formation in a glass melt depends on 
the intensity of the conversion of bonded silicon dioxide 
from the surface of dissolving quartz sand grains into the sili¬ 
cate melt and is estimated by the diffusion mass transfer co¬ 
efficient or, as it is usually termed, the coefficient of Si0 2 
diffusion. 

The values of Si0 2 diffusion coefficients have been ex¬ 
perimentally determined only for some compositions of so¬ 
dium-silicate glasses in limited temperature intervals that do 
not cover the regimes of melting impure sodium disilicate in 
industrial conditions. The experimental methods for deter¬ 
mining diffusion coefficients, as a rule, are technically com¬ 
plicated and labor-consuming, especially the methods based 
on direct measurement of thickness of the boundary diffusion 
layers formed around dissolving silica grains (the interfer¬ 
ence method, the radioactive isotope measurement method, 
x-ray analysis using an electron microprobe, etc.), since they 
involve special complicated and expensive measuring equip¬ 
ment [1-7]. 

At the same time, a significant spread is observed in dif¬ 
fusion coefficients of similar melt compositions. These dif¬ 
ferences can be accounted either by the difference of the 
methods of analysis used and their different accuracy, or by 
the degree of correctness of subsequent experimental data 
processing. The isolated experimental data available in the 
reference literature [8] are not fit for systematizing and ge¬ 
neralizing. 

In this context we proposed a method for the calculation 
of Si0 2 diffusion coefficients using kinetic curves of silica 
grain dissolution. 

Kinetic curves of grain dissolution in the melts of the 
Na 2 0 - SiO, system can be uniquely determined based on a 
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generalized dependence [9, 10]: 
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where m is the part of Si0 2 that has reacted; x is the relative 
duration of melting, sec; x c is the duration of complete melt¬ 
ing of the batch, sec; n is the power index, which is a func¬ 
tion of the initial diameter d 0 of Si0 2 particles in the batch, 
mm, and the melting temperature t, °C; and M is the silica 
modulus of the batch. 

Under the conditions of melting sodium-silicate glasses 
similar to industrial conditions, the silicate formation stage 
typically is completed already at x = 0.1; by that moment 
90 - 95% of initial SiO., quantity becomes chemical fixed. 

Upon passing to the glass formation stage, the rate of dis¬ 
solution of the residual grains sharply decelerates, stabilizes, 
and is determined only by diffusion mass transfer. 

It is known that the diffusion material flow (kg/sec) for 
the dissolution of a spherical particle in stationary conditions 
is found from the expression 

I=2ndDAy SlQ2 , (2) 

where d is the particle diameter, m; D is the molecular diffu¬ 
sion coefficient, m 2 /sec; Ay sic , 2 =Ys io , -y“ i02 is the con¬ 
centration head constituting the difference between Si0 2 
concentration on the particle surface and in the surrounding 
melt, kg/m 3 . 
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An ultimate saturation concentration of the melt for a 
given temperature y g i0 , is found from the solubility diagram 

of the system. 

The material flow is determined based on the material 
balance of a particle dissolving in the conditions of diffusion 
mass transfer and is calculated from the expression 
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Let us determine the time needed to complete the disso¬ 
lution of the residual sand grains: 

x est = 10,140 x (1 - 0.2) = 8112 sec. 

At the estimated point the homogenized composition of 
the silicate melt formed upon the dissolution of 96.3% initial 
silica contains 76.73% Si0 2 , which corresponds to M= 3.5. 

The Si0 2 concentration in this melt will be 


where p p is the density of the material of the particle, kg/m 3 . 

From joint analysis of expressions (2) and (3) we obtain 
a calculation formula for the diffusion coefficient: 
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TsiO, = 0.01(SiO 2 )p m = 0.01 x 76.73 x 2195 = 1684 kg/m 3 . 

The saturation concentration at a temperature of 1250°C is 
fsiOj = 0.01 (Si0 2 )p sm = 0.01 x 82.4 x 2193 = 1807 kg/m 3 . 


The following assumptions are made in the proposed 
method for calculating the diffusion coefficient: 

- the dissolving particles of quartz sand have a spherical 
shape; 

- the dissolution of particles of the same fraction in 
glass melt proceeds simultaneously according to the same ki¬ 
netic curve; 

- the delay of the process of homogenization of Si0 2 
content dissolved in the melt is not taken into account. 

Let us consider a particular example of calculating the 
diffusion coefficient using the generalized kinetic depen¬ 
dence. 

Let us take the reference data from an experimental melt¬ 
ing [11, 12], which investigated the process of batch melting 
of sodium silicate glass with M= 3.535 and the initial sand 
grain size 0.55 mm at a temperature of 1250°C. 

The generalized kinetic dependence (1) for the particular 
melting regime has the following form: 


The concentration head, accordingly, is equal to 
Ay si02 = 1807 - 1684 = 123 kg/m 3 . 


The values p m and p s m are determined using the Cohe- 
nen [8] data as values inverse to the specific volumes of the 
melt under the conditions considered. 

Using formula (4) we calculated the value of Si0 2 diffu¬ 
sion coefficient for the sodium silicate melt with M = 3.5: 


D = 


(0.183 xl0~ 3 ) 2 x 2650 
12x123x8112 


= 7.42 x 10 m 2 /sec. 


The diffusion coefficient value calculated by us agrees 
with known experimental data. In particular, the authors of 
[5], studying the dissolution of vitreous silica in a sodium si¬ 
licate melt with M= 2.0 at 1250°C, obtained the value 
D = 7.55 x 10 _ 12 m 2 /sec. 


m = 


-0.0232. 


_ _ X 

T_ 10,140’ 


REFERENCES 


since 

n = 0.05exp (3 x 0.055 2 ) x 1.25 “ 7 5 = 0.0232; 
x c = 6355 x 1.25 “ 10 83 ^ 915 M 2 - 717 = 10,140 sec. 

Let us arbitrarily fix a point on the kinetic cure segment 
x>0.1 (for instance x = 0.2), at which the dissolution pro¬ 
cess depends only on diffusion mass transfer, and then find 
the part of Si0 2 that has dissolved: 

m = o.2 0 0232 = 0.963. 

By that moment the residual sand grains will have the 
following estimated size: 

d est = ^(l-0.963)x0.55 3 =0.183 mm. 
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